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The  Hamiltonian  Stricture  of  the  Maxwe  11  -Vlasov  Equations 


Jerrcld  E.  Marsden*  and  Alan  Weinstein** 

Department  of  Mathematics 
University  of  California 
Berkeley,  CA  94720 

Abstract 

Morrison  [1980]  has  shown  that  the  Maxwell-Vlasov  and 
Vlasov-Poisson  equations  for  a  collisionless  plasma  can  be  written 
in  hamiltonian  form  relative  to  a  certain  Poisson  bracket.  We 
derive  another  (very  similar)  Poisson  structure  for  these  equations 
by  using  general  methods  of  sy^lectic  geometry.  The  main  ingredients 
in  our  construction  are  the  symplectic  structure  on  the  co-adjoint 
orbits  for  the  group  of  canonical  transformations  and  the  symplectic 

l 

structure  for  the  phase  space  of  the  electromagnetic  field  regarded 
as  a  gauge  theory.  Our  construction  shows  where  canonical  variables 
can  be  found  and  can  be  extended  to  Yang-Mills  gauge  theories. 


1.  Introduction 


In  this  paper  we  show  how  to  construct  a  Poisson  structure  for 
the  Maxwell-Vlasov  and  Vlasov-Poisson  equations  for  collisionless 


compare  our  structure  to  that  obtained  by  Morrison  [1980] 


#Tr  Sj!ft 

plasmas  by  using  general  methods  of  symplectic  geometry.  We  shall  /  ^  ^  0 

■ 
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We  consider  a  plasma  consisting  of  particles  with  charge  e  and 
■ass  ■  Moving  in  Euclidean  space  R3  with  positions  x  and 
velocities  v.  For  simplicity  we  consider  only  one  species  of 
particle;  the  general  case  is  similar.  Let  f(x,v,t)  be  the  plasma 
density  at  time  t,  and  E(x,t)  and  B(x,t)  be  the  electric  and 
magnetic  field.  The  Maxwell -Vlasov  equations  are: 


|i .  *.|£ .  .  ^).|f  •  o 

3t  3x  m'  c  /  3v 

(1.1) 

—  ~  -  curl  E 

c  3t 

A 

1 

i 

\ 

> 

(1.2) 

i||.curlB-|J  vf(x,v,t)dv 

I 

div  E  »  pf,  where  pf  •  e  j  f(x,v,t)dv 

f 

\ 

(1.3) 

div  B  »  0 

1 

Letting  c  -►  <»  leads  to  the  Vlasov-Poisson  equation: 

3t  3x  m  3x  3v 

(1.4) 

where 

"  Pf. 

(1.5) 

In  what  follows  we  shall  set  c  ■  1  in  the  Maxwell- Vlasov  system. 

The  Hamiltonian  for  the  Maxwell -Vlasov  system  is 

H(f,E,B)  ■  f  |mv2f(x,v,t)dxdv  ♦  f  i[E(x,t)2  ♦  B(x,t)2]dx  (1.6) 
Jx,v  L 


while  that  for  the  Vlasov-Poisson  equation  is 

H(f)  ■  [  inv2f(x,v,t)dxdv  ♦  j  [  <|>f  (x)pf  (x)dx  (1.7) 

'X,V  'x 

The  Poisson  bracket  used  by  Morrison  is  defined  on  functions 
F(f,E,B)  of  the  fields  f  ,E  ,B  by 


{F,G>  *  lx  vf{^^}dX<1V  +  lx  **  "  ^*CUrl  S  dX 

+  \  irSv  H  ■  frU  S?  dxdv 


lxv  fHMff-fl^ffdxdv 


(1.8) 


where  in  the  first  ten  {  ,  }  denotes  the  usual  Poisson  bracket  for 
functions  of  (x,v) ,  and  where  the  functional  derivatives  are  defined 
in  tens  of  the  usual  (Frichet)  derivative  by 


■  fi>v& 


■f'dxdv,  etc. 


For  the  Vlasov-Poisson  equation,  one  keeps  only  the  first  ten 
of  (1.8).  The  equations  (1.1)  and  (1.2)  (or  (1.4))  are  then 
equivalent  to 


F  -  {F,H} 


(1.9) 


with  H  given  by  (1.6)  (or  (1.7)  for  the  Vlasov-Poisson  equation  . 


Our  purpose  is  to  show  how  another  Poisson  structure  can  be 
constructed  by  a  general  procedure  involving  reduction  (Marsden 
and  Weinstein  [1974])  and  coupling  of  Hamiltonian  systems  to  gauge 
fields  (Weinstein  [1978b]).  Our  Poisson  bracket  differs  from  (1.8) 
in  that  the  last  integral  is  replaced  by 


f  6F  3f  r3_/5G\ 
I  fif  3v*L3v\6?' 

X  f  V 


B  dxdv 


<5G  3f  T3_/6G  \ 
6?  3v  Uvlfi?/' 


x 


B i dxdv 


(1. 


Both  structures  satisfy  the  Jacobi  identity  and  yield  the  correct 
equations  of  motion  for  the  hamiltonians  we  have  specified;  however, 
ours  is  constructed  by  general  principles  rather  than  being  ad  hoc, 
and  it  is  clear  where  the  canonical  variables  are  to  be  found.  In 
addition,  the  equations  div  E  *  p_-  and  div  B  *  0  arise  naturally  from 
the  gauge  symmetiy  of  the  problem  and  need  not  be  postulated  separately. 
We  believe  that  the  presence  of  the  two  symplectic  structures  might 
imply  the  existence  of  some  new  conserved  quantities.  11  Our  Poisson 
structure  fits  into  a  pattern,  special  cases  of  which  have  been 
found  by  others.  For  exaople,  Arnold  [1966]  showed  that  the  Euler 
equations  for  a  perfect  incompressible  fluid  are  a  Hamiltonian  system 
in  the  canonical  Poisson  structure  associated  with  the  group  of 
volume  preserving  diffeomorphisms  of  a  region  in  R3.  Using  Arnold's 
methods,  one  can  also  see  that  the  compressible  equations  are 
associated  to  the  semidirect  product  of  the  group  of  diffeomorphisms 
and  the  (additive  group  of)  densities  on  R3.  It  is  easy  to  check 


that  this  approach  yields  the  same  Poisson  structure  found  for 
perfect  fluids  by  Morrison  and  Greene  [1980]  and  ought  to  be  extendible 


s 


to  the  MiD  equations  by  the  methods  of  this  paper.  The  KdV  equation 
is  associated  with  the  Lie  algebra  of  the  group  of  canonical  trans¬ 
formations  in  the  work  of,  for  example,  Adler  [1979],  (We  recall 
that  there  is  a  standard  link  between  the  Maxwell-Vlasov  equations 
and  the  KdV  equation  --  see  Davidson  [1972]).  In  Ebin  and  Marsden 
[1970],  the  functional  analytic  machinery  required  to  fully  justify 
Arnold's  approach  was  given.  Itwas  proved,  for  example,  that  the 
volume  preserving  diffeomorphisms  form  a  true  infinite  dimensional 
manifold  which  is,  in  an  appropriate  sense,  a  Lie  group.  It  was  also 
shown  that  the  group  of  canonical  transformations  has  similar 
features,  but  no  physical  interpretation  was  given.  The  Vlasov 
equation  provides  one. 

All  of  these  clues  suggest  that  it  is  fruitful  to  find  a  more 
geometric  and  group-theoretic  framework  for  the  basic  equations  of 
plasma  physics.  Such  a  framework  is  provided  here.  We 
shall  not  deal  with  the  delicate  functional  analytic  issues  needed 
to  make  precise  all  the  infinite  dimensional  geometry,  nor  shall  we 
deal  with  questions  of  existence  and  uniqueness  (cf.  Batt  [1977]  and 
Horst  [1980]). 
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.  The  Poisson  Structure  for  the  Density  Variables 


We  begin  by  explaining  the  geometric  meaning  of  the  first  term 
in  (1.8),  |  +  |^]dxdv.  In  the  following  sections,  we  shall 

explain  the  term  for  Maxwell's  equations  (the  second  integral  in 
(1.8)),  and  then  finally  the  coupling  terms  (the  remaining  two 
integrals) . 

In  the  absence  of  a  magnetic  field,  by  normalizing  mass,  we 

can  identify  velocity  with  momentum;  hence  we  let  R6  denote  the 

usual  position-momentum  phase  space  with  coordinates 

(x1,x2,x3,p1,p2,p3)  and  symplectic  structure  Edxi  ~  dp^  (See 

Abraham  and  Marsden  [1978]  or  Arnold  [1978].)  Let  S  denote  the 

group  of  canonical  transformations  of  R6  which  have  polynomial 

growth  at  infinity  in  the  momentum  directions.  The  Lie  algebra  4 

of  S  consists  of  the  hamiltonian  vector  fields  on  R6  with 

polynomial  growth  in  the  momentum  directions.  We  can  identify 

* 

elements  of  4  with  their  generating  functions,  so  that 
4  •  d° OR6) ;  the  (right)  Lie  algebra  structure  is  given  by  the  negative 
of  the  usual  Poisson  bracket  on  phase  space.  (This  follows  from 
Exercise  4.1G  and  Corollary  3.3.18  of  Abraham  and  Marsden  [1978]). 

The  dual  space  4*  can  be  identified  with  the  distribution 
densities  on  R6  which  are  rapidly  decreasing  in  the  momentum 
directions;  the  pairing  between  h  €  4  and  f  G  4*  is  given  by 


The  generating  function  of  a  hamiltonian  vector  field  is  determined 
only  up  to  an  additive  constant,  pie  "correct"  group  S  is  really 
the  group  of  transformations  of  R°  x  R  preserving  the  1-form 
Epjdg.  ♦  dt  (Van  Hove  [1951]),  but  we  can  ignore  this  technical  point 
here  without  encountering  any  essential  difficulties. 
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<h,f)  ■  hf  dxdp  . 

(The  "density"  is  really  f  dxdp,  but  we  denote  it  simply  by  f.)  Now 
as  for  any  Lie  algebra,  the  dual  space  4*  carries  a  natural  Poisson 
structure  which  is  non-degenerate  on  the  co-adjoint  orbits  (see, 
e.g.  Guillemin  and  Sternberg  [1980]).  In  our  case  the  orbit 
through  f  €  4*  is 


0f  «{n*f  |  n  ^  4} 

and  so  the  Kirillov-Kostant  symplectic  form  on  0 ^  is  given  by 
o)f({f,h),{f,k})  «  <f,{h,k}> 

where  {f,h}  »  Ly  f  (Lie  differentiation)  represents  a  typical 

tangent  vector  to  0f  at  f.  (See  p.  303  of  Abraham  and  Marsden 

[1978]  —  two  minus  signs  have  cancelled  here.)  The  hamiltonian 

vector  field  Xp  on  0^  determined  by  a  smooth  function 
* 

F:  4  -*>  R  satisfies 

cof(Xp(f),{f,k»  -  dF(f)  »{f,k) 
for  all  k  6  4*.  We  claim  that 


(2.1) 


(2.2) 


(2.3) 


(2.4) 
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Indeed,  by  (2.3)  and  (2.2)  we  need  only  check  that 

-  <  =  dF(f) *{f,k}  =  <|£,{f,k}> 

But  (2.5)  holds  by  integration  by  parts.  In  fact,  the  following 
identity  is  of  general  utility: 

-  <f,{h,k}>-<{f,k},h> 


Thus,  the  Poisson  bracket  on  6*  is  given  by 


(F,G}(f)  »  wf(Xp(f),XG(f)) 


wf 


/  c  /6F  SGn  \ 

'  W* 


(by  definitic- 

(by  (2.4)) 

(by  2.2) 


We  have  proved: 

2 . 1  Proposition.  The  natural  Poisson  structure  on  the  dual  of  the 
Lie  algebra  of  the  group  of  canonical  transformations  is  given 


(F,G) 


f(g,  ffldxdp 


(2.5) 


(2.6) 


) 


(2.7) 
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Remarks .  1.  Notice  that  (2.7)  coincides  with  the  first  term  for 

the  Poisson  structure  (1.8)  if  p  is  replaced  by  y. 

2.  The  bracket  (2.7)  automatically  satisfies  the  Jacobi 
identity  since  it  coincides  with  the  Poisson  bracket  on  each  of 
the  symplectic  manifolds  0^. 

3.  (a)  If  f  is  a  delta  density,  0 ^  "coincides"  with 

(In  a  similar  way,  every  symplectic  manifold  is  a  co-adjoint  orbit.) 
For  f  a  sum  of  n  delta  functions,  0 ^  is  the  phase  space  for 
n  particles.  For  continuous  plasmas,  f  is  taken  to  be  a 
continuous  density,  in  which  case  0^  can  be  shown  to  be  a  smooth 
infinite  dimensional  manfiold. 

3.  (b)  If  f  is  a  density  concentrated  along  a  curve,  then 

is  identifiable  with  all  curves  having  a  fixed  action  integral. 

This  is  a  reduced  form  of  the  loop  space,  a  symplectic  manifold  used 
in  the  variational  principle  of  Weinstein  [1978a].  If  f  is 
concentrated  on  a  Lagrangian  torus,  then  0 ^  consists  of  lagrangian 
tori  with  fixed  action  integrals.  This  is  related  to  a  variational 
principle  of  Percival  [1979]. 

4.  By  using  an  appropriate  Darboux  theorem,  (see  Marsden  [1981], 
lecture  1),  one  can  show  that  0^  admits  canonically  conjugate 
coordinates . 

5.  The  Vlasov-Poisson  equation  is  a  hamiltonian  system  on  4* 
with  energy  function  given  by  (1.7).  If  f  evolves  according  to 
(1.4)  then  (1.9)  is  true.  This  is  a  direct  calculation,  already 
noted  by  Morrison  [1980].  More  can  be  learned  from  our  derivation 
of  the  Poisson  structure:  the  equation  (1.4)  is  tangent  to  each 


mmm 
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orbit  0 so  it  defines  a  hamiltonian  system  on  each  orbit.  This 
can  be  seen  directly  by  noting  that  (1.4)  can  be  written  in  terms  of 
ordinary  Poisson  brackets  as 

II  -{f,  «(f))  (1.4)  * 


where 


H(f)  -  ^mv2  +  ^4>f(x). 

Thus,  the  evolution  of  f  can  be  described  by 


where  f^  is  the  initial  value  of  f,  ft  is  its  value  at  time  t 
and  nt  e  S.  In  particular,  if  F  is  a  function  of  a  single  real 
variable  we  get  the  well-known  conservation  laws 

,  F(f.)  *  constant  in  time 
Jr6  t 

by  the  change  of  variables  formula  and  the  fact  that  each  n  €  S  is 
volume  preserving.  (These  conservation  laws  are  useful  in  proving 
existence  and  uniqueness  theorems  since,  as  in  the  case  of  two- 
dimensional  ideal  incompressible  flow,  they  lead  to  a  priori  L^- 
estimates.) 
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6.  In  Ebin-Marsden  [1970]  the  convective  term  v*7v  in  fluid 
mechanics  led  to  a  crucial  difference  between  working  spatially 
(in  the  Lie  algebra--the  "Euler"  picture)  or  materially  (on  the 
group... the  "Lagrange"  picture).  Here  there  is  no  such  term,  since 
it  would  be  given  by  (f,f),  which  vanishes. 

6.  Analogies  with  fluid  mechanics  raise  several  interesting 

analytical  problems:  (A)  if  6  times  a  dissipation  term  associated 
with  collisons  is  added,  do  the  solutions  converge  to  those  of  the 
Vlasov  equation  as  6  +  0?  (Analogous  to  the  limit  of  zero  viscosity). 
Standard  techniques  (Ebin-Marsden  [1970],  Kato  [1975])  can  probably 
be  used  to  answer  this  affimatively  for  short  time. 

(B)  Can  the  hamiltonian  structure  be  used  to  study  cha  otic  or 
turbulent  dynamics,  as  was  done  in,  for  exanple.  Holmes  and  Marsden 
[1981]? 

(C)  Is  the  time-t  map  for  the  Vlasov-Poisson  or  Maxwell-Vlasov 
equations  smooth?  See  Ratiu  [197V]  foT  a  discussion  of  why  this 
question  is  of  interest  for  the  KdV  equation. 
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§3.  Maxwell’s  Equations  and  Reduction 

Before  coupling  the  Vlasov  equation  to  the  electromagnetic 
field  equations,  we  shall  review  the  hamiltonian  description  of 
Maxwell's  equations.  The  appropriate  Poisson  bracket  for  the  E 
and  B  fields  (the  second  term  in  (1.8))  will  be  constructed  by  reduction 
(Marsden  and  Weinstein  [1974]). 

Let  P  •»  M  be  a  given  principal  S*  bundle  over  a  manifold  M. 

3  13 

In  our  case  M  *  R  and  P  ■  S  xS  .  Let  OL  denote  the  (affine) 
space  of  all  connections  for  this  bundle.  Let  G  denote  the  group 
of  gauge  transformations;  i.e.  bundle  automoxphisms  of  P  -►  M. 

i* 

Elements  of  G  may  be  denoted  e  where  M  -*■  F.  There  is  a 
natural  action  of  G  on  OL  given  by 

(e**,A)  M-  A  «•  d*  (3. 1) 


Consider  T *Ql  ,  the  cotangent  bundle  of  Cl  with  the  canonical 
synplectic  structure.  Elements  of  may  be  identified  with  pairs 

(A,E)  where  A  is  a  connection  and  E  is  a  vector  field  density 
on  M.  The  canonical  symplectic  structure  u  on  T*(X  is  given  by 


wKAj.EjMAj.Ej)) 


{h<va, 


-  Ej*A2)dx 


with  associated  Poisson  bracket  given  by 


(3.2) 


{F,G} 


f  6F  5G  5E  5G  \  . 

mtk  dx 
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Maxwell's  equations  in  tens  of  E  and  A  are  Hamilton's  equations 
on  T relative  to  (3.2)  for  the  haailtonian 


H(A,E) 


I E 1 2dx  ♦  j||dA|2dx 


(3 


Now  G  acts  on  A  by  (3.1)  and  hence  on  T *CL  .  It  has  a 
momentum  map  J:  TVC+q*  where  Cjj ,  the  Lie  algebra  of  G t  is 
identified  with  the  real  valued  functions  on  M.  The  momentum  map 
may  be  determined  by  a  standard  formula  (Abraham  and  Marsden  [l978. 
Corollary  4.2. ll])  to  be: 


J(A,E)*$ 


|p*d<fc>d*s-j;div 


A\ 


(3 


Thus  J_1(0)  ■  {(A,E)  €  T*  |div  E  ■  0)  .  By  a  general  theorem  on 
reduction  (Marsden  and  Weinstein  [1974]),  the  manifold  J_1(0 )/G 
has  a  naturally  induced  symplectic  structure. 


3.1.  Proposition.  We  may  identify  J-1(0)/G  with 
M  «  {(E,B)  |  div  E  ■  0,  div  B  ■  0).  The  Poisson  bracket  on  M  is 
given  by 


{F,G} 


6g  6G 
3b  "  3e 


curl 


(3 


Maxwell’s  equations  in  the  vacuum  are  Hamilton’s  equations  far. . 


H(E,B) 


j  J |E {2dx  ♦ 


(3 


Proof.  The  identification  between  J  (0 )/G  and  M  is  by  the 


equation  B  -  curl  A.  (We  assume  M  has  trivial  first  cohonology, 
as  is  the  case  here;  for  Yang  Mills  fields,  B  is  the  "curvature" 
of  A.)  Let  F,G:  M  -*•  F.  We  aay  regard  F  and  G  as  G- invariant 

a  a 

functions  F,  G  on  T*.Y  by  B  ■  curl  A.  Since  Poisson  brackets  push 
down  under  reduction  we  have 

{F,G)0  ■  {F,G}  C 

where  0  aeans  the  induced  function  on  M.  Now  in  the  canonical 
structure  on  T<t  we  have 


.  ^  a, 

{F.G} 


/A  A  A  A 

f6F  £G  6G  6F  . 
>31  TB  '  IA  / 


dx 


( 


The  chain  rule  and  the  definition  of  functional  derivatives  give  the 
identity 


jyjj»A*dx»  curl  A  'dx-  -  jA’curl  jg-dx  ( 

Substitution  of  (3.9)  in  (3.8)  gives  (3.5).  The  rest  is  readily 
checked.  ■ 

This  foraalisa  generalizes  readily  to  Yang-Mills  fields  and 
to  these  fields  coupled  to  gravity;  see  Arms  [1979]. 
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§4.  A  General  Construction  for  Interacting  Systems 

The  work  of  Weinstein  [1978b]  on  the  equations  of  motion  for  a 
particle  in  a  Yang-Mills  field  uses  the  following  general  set-up. 

Let  it:  8  -*•  M  be  a  principal  G-bundle  and  Q  a  haailtonian  G-space 
(or  a  Poisson  manifold  which  is  a  union  of  hamiltonian  G-spaces). 

Then  G  acts  on  T*8  and  on  Q,  so  it  acts  on  Q  *  T*8  .  This  action 

has  a  momentum  map  J  and  so  may  be  reduced  at  0: 

(Q  *  T*8)0  -  J”1 (0) /G 

The  reduced  manifold  carries  a  syntactic  (or  Poisson,  if  Q  was 
a  Poisson  manifold)  structure  naturally  induced  from  those  of  Q 
and  T*8. 

For  a  particle  in  a  Yang-Mills  field  one  chooses  8  to  be  a 

G-bundle  over  3-space  M  and  Q  a  co-adjoint  orbit  for  G  (the 

internal  variables).  The  hamiltonian  is  constructed  using  a 
connection  (i.e.  a  Yang-Mills  field)  for  8.  For  electromagnetism 
G  ■  S1  and  Q  *  (e)  is  a  point. 

For  the  Vlasov-Maxwell  system  we  choose  our  gauge  bundle  to  be 

8  -  *  -*•  M 

where  M  ■  { (E,B)  |  div  B  ■  0}  ,  with  G  »  G  the  gauge  group 
described  in  the  previous  section. 

We  now  choose  S  to  be  the  group  of  canonical  trans format ions 
of  T*M  (»F6).  We  can  let  Q  be  eitner  the  symplectic  manifold 


T*S  or  the  Poisson  manifold  6*.  It  remains  to  specify  an  action 
of  G  on  T*S  .  We  set  e  ■  1  and  let  G  act  on  S  by 
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($,n)  **■  n  ♦  d^oiron  (4.i) 

where  it:  F^  F"*  is  the  projection  (x,p)  »*•  x;  i.e.  we  translate  the 
momentum  space  at  x  by  d$(x).  A  simple  computation  gives: 


4.1.  Lemma.  The  action  of  G  induces  an  action  on  4*  given  by 


(e^.f)*-*  for 

where  €  S  is  translation  by 
for  this  action  is  given  by 


d$ 

d*  . 


The  momentum  map  J:  4*  G* 


•I 


f(x,p)<J>(x)dxdp 


To  construct  the  interaction  space  we  need  to  compute  the 
momentum  map  for  the  corresponding  action  of  G  on  4*  x  T*dt  given 
by 


(el*(f.A.E))  m.  (for^.A  ♦  d*,E) 


(4.2) 
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J(f,A,E)-4> 


f 

l 


f(x,pH(x)dxdp  ♦  E»d$dx 
x.P  J 


Pf(x)<Kx)dx  -  div  E  ♦(x)dx 

X  'X 


Again  this  is  a  straightforward  calculation.  Our  reduced 
(interaction)  space  J'*(0)/G  may  now  be  identified  with 

I/M  -  {(f,B,E)  |  div  E  *  pf,  div  B  -  0) 


where 


(4. 


J_1(0)  »  { (f,A,E)  |  div  E  -  pf}. 
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§5.  Computation  of  the  Poisson  Structure 

By  reduction,  the  space  VM  defined  by  (4.3)  has  a  natural 
Poisson  structure;  we  now  compute  it. 


5.1.  Theorem.  The  Poisson  structure  on  VM  is  given  by 


{F,G} 


*  f„(H curi  w  •  H curi  S)dx 

*  [|  f(lp  frXfe  4’1<iiv  Hr)d’‘d>’ 

L  X»P 

’  (  £(lp  l?fe(A’1<iiv  H)dxdp] 

*»p 


Proof.  Having  identified  J 
that  the  natural  projection  ir:  J-1 


_1C0)/6 

(0)  VM 


with  VM  ,  one  checks 
is  given  by 


ff(f,A,E) 


curl  A,  E) 


where 


d>A  ■  -A"1(div  A) 


and,  as  before,  t ..  is  translation  by  d$. 

A 


(5.1) 


(5.2) 


(5.3) 
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Next,  define  a  projection  map 

X:  -4*  x  T *0T.  -►  J_1(0)  ;  X(f,A,E)  •  (f ,A,E  ♦  7$f  £)  (5.4) 

where  £  *  A-1(p£  -  div  E) . 

Now  if  F,G  :l/M  ■*>  F,  we  can  extend  F,G  to  maps  on  4*  *  7*01. 
by 


F  -  FoIIoX,  G  -  GoIIoX  (5.5) 


If  we  denote  the  map  ll<>x  by 

(f»A,E)  (?,§,g) 


we  have,  from  (5.2)  and  (5.4), 


B 

E 


m  fOT 

d^A 

■  curl  A 


■  E  ♦  7$ 


f.E 


(5.6) 


Since  J-1(0)  is  co- isotropic,  Poisson  brackets  can  be  computed 
in  4*  *  T*X,  and  the  answer  is  independent  of  how  F  and  G  are 
extended  to  4*  x  T *CL  .  We  chose  the  specific  map  X  to  effect 
the  extension  explicitly. 
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Using  the  definition  of  functional  derivatives  and  the  chain 
rule,  one  finds  the  following  relationships 


S  -  -  A-1  div  " 

51  5?  5E 


<SF  f  3  A-l..  / 6F  3f\j 

SA--J 


.  5F 

curl  — 
6B 


(5.7) 


6F  6F  3  <SF 

Jf  '  T7  ‘  5T  Ax  dlv-  - 


5E 


5fi 


Substituting  (5.7)  into  the  Poisson  structure  on  e*  *  T*X  obtained 
from  that  on  4*  (see  2.7)  plus  that  on  T*2l  (see  3.2),  i.e.  into 

{F,G}  •  |  f{||,  ||}dx<*P  ♦  |  i.fjff  -  Msf  dx  (5-8) 

'X,p  1  X  ' 

we  obtain  formula  (5.1).  ■ 


Note:  In  carrying  out  the  substitution  of  (5.7)  into  (5.8),  some 

3  -1 

cancellation  occurs  since  ^  &x  div^V  is  the  gradient  part  of  V, 

2  3 

which  is  L  -orthogonal  on  K  to  divergence- free  vector  fields. 

The  Poisson  structure  (S.l)  automatically  satisfies  the  Jacobi 
identity,  since  it  is  the  structure  associated  with  reduced  manifolds, 
which  are  syqplectic. 


The  haailtonian  for  the  Vlasov-Maxwell  equations  in  the  momentum 
space  representation  is  the  function  H  on  6*  x  T *CL  given  by 

H(f,A,E)  «  f  ^(p-A)2f(x,p)dxdp  ♦  if  (E2  ♦  (curl  A)2)dx  (5.9) 

Jx,p  Jx 

This  haailtonian  is  gauge- invariant;  i.e.  it  is  invariant  under  the 
action  of  G  given  by  (4.2).  Thus  H  produces  a  well-defined 
haailtonian  on  l/M  .  We  find  that 

2 

H0(f,B,E)  «  f  E-  f(x,p-V  x  A_1B)dxdp  +  I  f (E 2  ♦  B2)dx  (5.10) 

Jx,p  ' 

5.2.  Theorem.  The  hamiltonian  (5.10)  in  the  Poisson  structure  (5.1) 
yields  the  Vlasov-Maxwell  equations  in  momentum  representation. 


This  can  be  checked  directly;  however,  to  facilitate  comparison 
with  Morrison's  results  we  shall  transform  our  Poisson  structure  to 
produce  a  velocity  space  verson  of  5.2. 
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§6.  Transformation  to  Velocity  Representation 

To  compare  our  bracket  (5.1)  with  Morrison's  (1.8),  we  transform 
(5.1)  to  velocity  space.  Let  us  denote  the  space  defined  by  (4.3) 
by  l/Mp  to  emphasize  that  f  is  a  density  on  position -momentum 
space.  We  write  l/My  for  the  corresponding  space  (4.3)  for  f  a 
density  on  position-velocity  space.  We  shall  take  e  «  1,  m  =  1 
for  simplicity.  Define  the  transformation 


P  :  >  VMp 


by 


(f,B,E)  w  (f,B,E) 


where 


and 


f(x,p)  «  f(x„p  +  V  *  A~*B(x)  ) 
B  *  B 
E  *  E 


(6.1) 


6.1.  Theorem.  The  transformation  (6.1)  transforms  the  bracket  (5.1) 
to  (1.8)  with  the  last  two  terms  replaced  by  (1.10).  Both  brackets 
yield  the  correct  equations  of  motion;  i.e.  (l.l)-(l. 3)  imply  that 


(1.9)  holds. 
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Sketch  of  Proof.  We  make  these  remarks:  (a)  Transforming  the 
first  term  of  (5.1)  produces  the  first  term  of  (1.8)  plus  the 
term  (1.10).  This  is  seen  using  the  chain  rule  and  the  equation 
p  ♦  7xA_1B  *  v. 

(b)  The  second  and  third  terms  of  (5.1)  both  produce  extra 
terms  via  the  chain  rule.  These  recombine  to  give  the  second  and 
third  terms  of  1.8’. 

(c)  In  the  special  case  in  which  G  is  the  Hamiltonian  (1.6), 
both  (1.10)  and  the  last  teim  of  (1.8)  give  the  same  expression.  This 
is  why  both  structures  give  the  same  equations  of  motion.  ■ 
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§7.  Additional  Remarks 

(A)  Poisson  structures  may  be  viewed  as  bundle  maps  taking 
covectors  to  vectors.  (In  this  guise,  they  are  called  "cosyaplectic 
structures.")  Viewed  this  way,  Morrison's  bracket  is  the  map 

(f* ,B* ,E*)  h.  (6f ,6B,6E) 


given  by 


Sf  -  -  ||*E*  -  ||*  (v  x  B*) 


6B 


fp(H‘  *  *) 
LS 


f*dv  -  curl  E* 


6E  =  f*dv  +  curl  B* 

P 


while  our  cosymplectic  structure  is  given  by 


(7.1) 


6f  *  -{f,f*  ♦  A'1  div  E*} 
6B  ■  -  curl  E 

«E  «  P  j  f*  ||dp  ♦  curl  B* 


(7.2) 


in  momentum  representation,  where  P  is  the  projection  to  the 
divergence  free  part.  In  velocity  representation,  (7.2)  becomes 
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Sf  •  -{f,f*  ♦  A'1  div  E*}  -  x  b) 

SB  =  -  curl  E 

5E  ■  |  f*dv  +  curl  B* 

Thus  (7.1)  and  (7.3)  differ  only  in  the  terms 

3f* 

v  x  B*  verses  x  B 
3v 


(7.3) 


From  (7.2)  we  see  explicitly  that  6f  is  tangent  to  co-adjoint 
orbits  in  S.  (This  fact  is  obscured  in  velocity  representation.) 

On  the  other  hand,  when  Morrison's  structure  (7.1)  is  transformed 
to  momentum  space,  6f  is  not  tangent  to  co-adjoint  orbits  of  S. 

(B)  A  "cold  plasma"  may  be  defined  as  one  for  which  f  is  a 

6  measure  supported  on  the  graph  of  a  vector  field  p  »  0(x).  This 
property  persists  as  f  evolves  by  composition  with  a  canonical 
transformation.  In  fact,  the  property  that  0  is  curl-free  is  also  main¬ 
tained,  since  this  corresponds  to  the  graph's  being  a  lagrangian 
submanifold.  After  a  long  time,  the  submanifold  may  no  longer  be 
a  graph.  This  is  the  "shock"  phenomenon,  leading  to  multiple 
streaming  (Davidson  [1972]  .)  We  remark  that  Maslov  ([1976],  p.  44) 
has  already  noticed  this  evolution  of  lagrangian  submanifolds  for  the 
Vlasov- Poisson  equation. 

(C)  We  would  like  to  understand  in  general  terms  the  contraction 
of  one  hamiltonian  system  to  another.  Exanqiles  are  the  passage  to  the 
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restricted  three  body  problem  from  the  full  three  body  problem  and 
the  limit  c  ■*  «  to  get  the  Vlasov-Poisson  equation.  It  would  also 
be  of  interest  to  realize  both  the  Vlasov- Maxwell  and  MID  equations 
as  limiting  cases  of  a  grand  hamiltonian  system  constructed  from 
the  Boltzmann  equation. 

(D)  We  have  remarked  that  our  formalism  readily  generalizes 
to  Yang-Mills  interactions.  Is  such  a  hamiltonian  structure  useful 
in  nuclear  physics  for  Yang-Mills  plasmas? 
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